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Abstract. - We present a novel approximation scheme to describe the influence of a harmonic 
bath on the dynamics of a two-level particle over almost the whole regime of temperatures and 
coupling to the environment, for a wide class of bath spectral densities. Starting from the exact 
path-integral solution for the two-level system density matrix, effective intra-blip correlations are 
fully included, while inter-blip and blip-sojourn interactions are considered up to first order. In the 
proper regimes, an excellent agreement with conventional perturbative approaches and ab-initio 
path-integral results is found. 
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Introduction. — The problem of a two-level system 
(TLS) suffering from environmental decohering effects is 
ubiquitous to many physical and chemical situations [1-3]. 
Standard examples involve electron and proton transfer 
reaction in condensed phases [4] , defect tunneling in met- 
als [5] or tunneling systems in glasses [6,7]. Recently sev- 
eral realizations of TLSs have been experimentally demon- 
strated in superconducting [8] and semiconducting [9] de- 
vices as possible unit (quantum bit) for future quantum 
computers. In these solid state systems, decoherence is a 
major obstacle towards the realization of a usable quan- 
tum computer [10-12]. Hence, a proper understanding of 
dissipation over a broad parameter regime is of outermost 
importance. 

' For a description of the dissipative dynamics the spin- 
boson model, in which the TLS is bilinearly coupled to a 
harmonic bath, is very frequently used. It reads [1-3] 



H(t) = -W)a z 



H 



B ■ 



(1) 



The basis states |R) and |L) are the localized eigenstates 
of the "position" operator <r z , A describes the coupling 
between the two-states due to tunneling, and e(t) is an ex- 
ternal control field. The Hamiltonian = J^i ^i{b\oi + 
1/2) represents a bath of bosons, and the collective vari- 
able X = c i(bi + Sj)/2 describes the bath polarization. 
Despite the huge amounts of works on the subject [1-3], 



the existing schemes for a portrayal of the time-evolution 
of the TLS reduced density matrix mostly reduce to two 
main roads of approximation. On the one hand the so 
termed noninteracting-blip approximation (NIB A) [1,2], 
or equivalent projection operator techniques [13] based 
on an expansion to leading order in the tunneling ma- 
trix element A, has been proved to be successful in the 
regimes of high temperatures and/or strong friction. On 
the other hand the weak coupling and low-temperature 
regime, where NIBA fails for an asymmetric TLS, is typ- 
ically tackled within an expansion to lowest order in the 
TLS-bath coupling. In this latter case path-integral meth- 
ods [14, 15] as well as the Bloch-Redfield formalism are 
used [16] (the two methods have been demonstrated to 
yield the same dynamics for weak Ohmic damping [17]), 
or a Born approximation [18]. To date, only numerical ab- 
initio calculations [12, 19-21] can provide a description of 
the TLS dynamics smoothly interpolating between a weak 
and a strong coupling situation. 

In this work, we present an interpolating approximation 
scheme, enabling to describe the weak and strong coupling 
regimes in a unique scheme. We call it weakly-interacting 
blip approximation (WIBA), within which the dynamics 
of the population difference (<x 2 )t = P(t) is 



P(t) 



dt'[K a (t,t')-W{t,t r ) + K s {t,t r )P(t')} . (2) 
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Fig. 1: Generic path with 2n = 6 transitions at flip times 
ti, ti, ..t2n- The system is in an off-diagonal state (blip) of 
the reduced density matrix in the time intervals Tj = tij — 
and in a diagonal state (sojourn) at times t2j+i — t2j- 
The interactions S^j^j-i, R2j,2j-i and Yjj-i (intra-dipole and 
blip-preceeding-sojourn interactions), Eq. (7), are symbolized 
by the wiggled lines (blue and magenta online, respectively). 
The double-dashed lines denote the inter-dipole interactions 
A.,^, while the bold-dotted lines are the remaining blip-sojourn 
interactions Xj±, cf. Eq. (ff). 

The irreducible kernels K^- s / a \ W entering this generalized 
master equation are neither perturbative in the tunneling 
matrix A nor in the TLS-bath coupling, and are given in 
analytical form in (18), (19) and (20) below. By compar- 
ing the predictions of the WIBA with known perturbative 
results as well as with exact ab-initio calculations, we show 
that the WIBA well describes the TLS dynamics over the 
whole regime of temperature and environmental coupling. 

In the spin-boson model all the effects of the bath on 
the TLS are captured by the spectral density G{lo) = 
nhr 2 cf 8(10 — LOi). In the following we shall consider 
a class of spectral densities with a continuous spectrum: 



G{u) = 28 s uj 1 - s uj s e 



with S s being a dimensionless coupling parameter, w p h a 
characteristic phonon frequency, and to c the bath cut-off 
frequency. Thus, (3) encompasses the commonly consid- 
ered Ohmic spectrum (s = 1) [1-5,13,17-19,21,22], with 
Si = a being the so-called TLS Kondo parameter, and the 
super-Ohmic case [1,2,6,7,11]. The applicability of the 
WIBA to other classes, as e.g. structured baths [4,12], 
will be discussed elsewhere. 

Exact path-integral formulation. To start with, 
we assume a factorized initial condition at time t = with 
the particle having been held at the site \R) (a z = +1) 
from time to = — 00 till t = 0, and with the bath in thermal 
equilibrium. Then the exact formal solution for P(t) can 
be expressed in terms of a real time double path integral 
over forward ct(t) and backward ct'(t) spin paths [2, 3] 
with piecewise constant values ±1. Upon introducing the 
linear combinations r]{ T )/£,( T ) = [c(r)±cr'(r)]/2, one finds 



P(t) = J V£Vr)A[£,ri]exp{$[Z,r)}} , 



(4) 



where A is the path weight in the absence of the bath 
coupling. A generic double path can now be visualized 
as a single path over the four-states of the reduced den- 
sity matrix, characterized by {j]{r) = ±1,£(t) = 0) and 
(v( T ) = 0,£(t) = il)- The time intervals spent in a di- 
agonal (£(t) = 0) and off-diagonal (ry(r) = 0) state are 



A 2 ,i 



Fig. 2: Irreducible kernel K {s) {t,t') in the NIBA (a), in the 
extended-NlBA (blue online) (b) and in the WIBA (c). The 
single-dashed lines are the linearized blip-blip interactions be- 
tween the first and last dipole. The inner bubble denotes the 
infinite sum of extended-~NTBA diagrams yielding the symmet- 
ric part of P(t) within the extended-NIBA, denoted . 



dubbed "sojourns" and "blips", respectively [1]. Due to 
the initial condition, the path sum runs over all paths 
with boundary conditions £(0) = £(t) = and 77(0) = 1, 
r](t) = ±1. Environmental effects arc in the functional 

$[e,r?]= I dt 2 [ 2 dt 1 i(t 2 )\s 2 , 1 i(t 1 )+iR2,ifi(t 1 )l (5) 
Jo Jo 1 J 

with the bath correlation function Q = S + iR being 



Q(t) 



dto 



G(u) 



coth 



2k B T 



(1 — cos tot) + i sina;t 



(6) 

and Qjk : = Q{tj — tk)- For a generic path with 2n 
transitions at times tj, j = l,2,...,2n, one finds £(t) = 
E?=i ^ S ( T ~ tj) and fi(r) = £)f= %<5(t - tj). Here is 
770 = 1 due to the initial preparation and £j = ±1, r\j = ±1 
for j > 0. Because ^23 = ~£2j-ij the influence function 
in (5) becomes $00 = $%L, hps + *2L ( Fi S- Thc 
function 3>}™ tra bps describes intra-blip and blip-preceeding 
sojourn correlations, and reads 



<I> 



00 

intra, bps 



n 

■n 



S2j,2j-1 l-V,., , 
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(8) 
(9) 

(10) 



where we split Aj,j„i = R 2 j, 2 j-\ + with 

Xjj-l = #2j-l,2j-2 - R-2j,2j-2 ■ 

Moreover, the functional &{2ter accounts for inter-blip 
and blip-sojourns interactions [1,2] 

n j-1 n j-2 



$ 



(») 



E E && A *.* + 1 E E ti^XiM • ( n ) 



3=2 k=i 



3=2 fe=0 



The function Aj.k contains the blip-blip interactions be- 
tween the flip pairs {j, k}, while the blip-sojourn interac- 
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tion Xj k yields a phase factor. To be definite, for k > 0, 

Aj,fe= S2j,2k-1 + S2j-l,2k — Sij t 2k — S2j-l,2fe-lj (12a) 
Xj t k = R2j,2k+1 + #2j-l,2fc — Rlj,2k — Rzj-lflk+l ■ (12b) 

The correlations -Xj- o depend on the initial prepara- 
tion [2]. The summation over the path histories then re- 
duces to an expansion in the number of tunneling transi- 
tions yielding formally exact, but practically untractable, 
equations for P(t) of the form (2) [14]. 

Known and novel approximation schemes. To 

find appropriate approximation schemes to the TLS dy- 
namics, let us start from the familiar non-interacting-blip 
approximation (NIB A) [1,2]. Within the NIB A, one sets 
3>hrter = 0, namely the inter-blip correlations Aj.k and 
the blip-sojourn interactions Xj t k (k ^ j — 1) are ne- 
glected. The blip-preceeding-sojourn interactions Yjj-i 
in Eq. (10) are neglected as well. Hence, Xj.j_i reduces 
to Xj t j-i ~ R2j,2j-i- The influence function (7) then 
splits into individual influence factors depending only on 
the dipole length Tj :— t 2 j — t 2 j-i- The dynamics is thus 
described by (2) with NIBA kernels corresponding to the 
one-dipole irreducible contributions (Fig. 2a), 

i^(M') = A 2 C(t-i')cos[C(M')] , 
KS(t,t') = A 2 S(t-t') S m[at,t')], 

with C(M') = Jl dt"e(t"), intra-blip contributions C(t) = 
e- s ^cos[R{t% and S{t) = e - s ^ sm[R(t)}. Here, W N = 
0. The kernels are of lowest order in the tunneling matrix 
A but are non-perturbative in 6 S . Due to the simplicity 
of the kernels (13), the NIBA has been a very popular 
approximation so far. For sub-Ohmic damping, s < 1, 
NIBA is expected to be a valid approximation for all tem- 
peratures with the TLS exhibiting incoherent dynamics 
even for very small coupling 6 S . For Ohmic and super- 
Ohmic damping, NIBA is expected to be a good approx- 
imation only at high enough temperature and/or strong 
damping [2]. However, its limit of validity are not clearly 
defined. The NIBA is known to fail at low temperatures 
and weak coupling for an asymmetric TLS for Ohmic and 
super-Ohmic damping, because the dipole-dipole correla- 
tions Aj t k contribute already to terms which depend lin- 
early on the spectral density G(uj). For example, in the 
case of a TLS with static asymmetry e(t) = Eo, NIBA 
predicts the unphysical asymptotic limit 

p- = -tanh(^) ! (14) 

implying localization of the TLS (P§f = —1) at zero tem- 
perature even for infinitesimal asymmetries. In order to 
overcome the NIBA shortcomings, a weak-coupling ap- 
proximation (WCA) has been proposed in [2, 14, 15] with 
WCA kernels being linear in 6 S and nonperturbativc in A. 
Within the WCA, the TLS dynamics shows damped coher- 
ent oscillations with a renormalized energy splitting Ml, 



with Q 2 = A 2 ff [1 — 2Reu(iE/h)] + £q, towards the equi- 
librium value 

^WCA = -^tanh(^). (15) 

Here, the frequency shift is related to the frequency in- 
tegral u(z) ~ | J °° duj ^^l-i [coth(fc/?o;/2) — 1]. Moreover, 
E = hyj A 2 S + £q and the effective bath-renormalized tun- 
neling coupling A c ff for the cases s > 1 reads [2] 

A off = A[r(l - 2a) cos(Tra)] ^t^t (A/uj^^ , s = l, 
A cff = Acxp^r^-l)^/^)^ 1 ], s>l, (16) 

with r(z) the Gamma function. Finally, the relaxation 
r r and dephasing rates are given by the perturba- 
tive (in S s ) expressions r r = (-nh 2 A 2 S /2E 2 )A(E /h) and 
T = r r /2 + {irh 2 el/2E 2 )A(Q), where the spectral func- 
tion A(uj) = G(u>) coth(fiw/2fcBT) is related to emission 
and absorption of a single phonon. 

To smoothly bridge between the high and low T lim- 
its, let us now start to discuss a more refined approx- 
imation, which we call extended-NIBA (Fig. 2b). As 
in NIBA, $i ntor is neglected, while the approximation 
on the blip-preceeding sojourn interactions Xjj—i is im- 
proved, considering also Y, j-i in an effective way. Specif- 
ically, expanding Yjj-i = 2 + R23 \:>. 1 up 
to first order in the blip lengths Tj, we set Xjj-i « 
R(rj) — Tji?2j-i,2j-2- As a result, the extended-NIBA ker- 
nels K s J^(t,t') have the same form as (13) with C(t) — > 
C'(t) = e- s ^ cos[R{t)}. Here is R(t) = R(t) - tR(t). 
Moreover, W c n = K^n — K^. A comparison between 
NIBA and extended-NIBA, as well as with other approxi- 
mation schemes discussed below, is shown in Figs. 3a - 3d. 
In Figs. 3b - 3d we also show results obtained with the nu- 
merical ab-initio path-integral approach QUAPI [20] . The 
short-time dynamics is always well approximated by the 
NIBA (extended-NIBA). At long times, however, correla- 
tions neglected in the NIBA become relevant. In partic- 
ular, already at intermediate temperatures and damping 
(Figs. 3c, 3d) the extended-NIBA correctly reproduces 
the QUAPI results while NIBA fails. At low T and small 
damping, both NIBA and extended-NIBA fail to repro- 
duce the correct long time dynamics (Fig. 3a), as pre- 
dicted e.g. from the WCA. 

To bridge between the moderate damping situation well 
described by the extended-NIBA and the extremely under- 
damped case we observe that, for spectral densities of the 
form (3), the blip-blip interaction terms Aj t k as well as 
the blip-sojourn terms Xj t k (k ^ j — 1) arc intrinsically 
small compared to unity. Hence, we propose a novel ap- 
proximation scheme, which we call the weakly-interacting 
blip approximation (WIBA). Within the WIBA, the full 
^ intra bps * s retained as in the extended-NIBA and one ex- 
pands the influence functional exp {$& er } up to linear 
order in the blip-blip and blip-preceeding sojourns inter- 
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Fig. 3: Time evolution of the expectation value P(t) at low/moderate temperatures k^T < E for several values of the Ohmic 
coupling parameter a. Full lines depict the WIBA, dashed lines the NIBA, dotted lines the extended-'NIBA, the dot-dashed 
ones are results for the weak-coupling approach (WCA) while the lines with bullets are the ab-initio QUAPI predictions. All 
quantities are expressed in units of A. At low damping and temperatures, Fig. 3a, the TLS exhibits damped coherent oscillations 
towards the asymptotic value -Pwca- As the damping is increased, the oscillations are more strongly damped, see Figs. 3b, 3c, 
3d. In the chosen regime of parameters, the equilibrium value P°° is neither well described by P^ nor by Pwca- 



actions Aj^ and Xj^- In other terms, 

exp {$(»)} « exp {$^ ra , bps } (l + §«) . (17) 

Within an expansion in the number of tunneling transi- 
tions, the lowest order self-energy corresponds to the ex- 
tended-NIBA, while higher order terms describe a set of 
blips in which the first and last blip are interacting, Fig. 
2c. Summing up the higher contributions, the WIBA ker- 
nels, neither perturbative in A nor in 5 S , read 

#w(*4,*i) = -Krf*(*4,*l) 

- A 4 f 4 dt 3 f 3 dt 2 C'{t A - t 3 )sin[C(t4,i3)]i 5 e N (*3 - h) 
x A2.iC'{t 2 -t 1 ) sin[C(i 2 ,*i)] , (18) 



#w(*4,*l) = #Ssr(*4)*l) 

- A 4 / dt 3 / dt 2 C'{U - h) Sin[C(i4,t3)]-PeN(*3 - h) 

x [-A 3l i5(t a - h) + X 2 , C(t 2 - h)] cos[C(i 2 , h)} . (19) 
Moreover, 

- A 4 f * dt 3 f 3 dt 2 C'(t 4 - ts) sin[C(k, ta)]iSr(*s - *a) 

Jt x Jt! 

x [A 3l i<yC(* a - *i) - X 2 . S(t 2 - h)] sin[C(t 2 , ti)] > ( 20 ) 

where 8C = C — C. Moreover, P^{t) is the symmetric 
part (in £o) of P o n(0 within the extended-NIBA. Thus, 
at high temperatures, where the blip-blip interactions are 
negligible, the WIBA kernels reduce to the extended-NIBA 
ones. By expanding the WIBA kernels to first order in <5 S 
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Fig. 4: Time evolution of P(t) for super-Ohmic damping (coupling parameter S3). Here the ratio w c /w p h = 8 is kept constant 
in both panels and we set T = 0.1, so = 1 (in units of A). Full lines are the WIBA predictions, dot-dashed lines are results of 
the WCA, dashed lines represent the NIBA predictions, dotted lines denote the extended-NlBA dynamics and finally bulleted 
lines are results of QUAPI. 



and approximating Xj^—i to B,2j-2j-i, the weak damping 
kernels in [2, 14] are recovered. 

Ohmic damping. — As a benchmark for the WIBA, 
we consider the evolution of the population difference P(t) 
for the important case of Ohmic damping. In Fig. 3, the 
Ohmic case (oj c = 50A and e(t) = Eq = A) is shown. An 
excellent agreement is found for weak damping and tem- 
peratures (Fig. 3a) between WIBA and WCA, whereas 
the extended-NIBA matches the NIBA and predicts the 
wrong asymptotic limit P™. As the coupling is increased 
(or by raising the temperature), the WCA is expected to 
fail. However, as for the NIBA, the limits of validity of 
the WCA are not clearly defined. Indeed, Figs. 3b to 3d 
show an intermediate parameter regime where both ap- 
proximations fail, since dipole-dipole interactions as well 
as two-phonon processes are relevant. Comparison with 
results from QUAPI shows that the short time dynamics 
is well approximated by the NIBA (WIBA) . At intermedi- 
ate and long times, the WIBA reasonably well approaches 
QUAPI and its asymptotic value. From a comparison with 
QUAPI, we notice that the higher order dipolc correla- 
tions neglected in the WIBA yield a larger dephasing rate 
than predicted from WIBA. In particular, QUAPI pre- 
dicts a complete suppression of the coherent oscillations 
already at a = 0.25, T = 0.01. An interesting case is 
shown in Fig. 3b, with a small-to-intermediate value of 
the coupling strength (a = 0.1), where the extended-NIBA 
slightly moves from the NIBA towards the WIBA predic- 
tions, reaching an intermediate asymptotic value (see in- 
set). 

Super-Ohmic damping. — Let us now consider the 
predictions of the WIBA in the super-Ohmic case (s = 3). 
Since S(t ~ A -1 ) differs only little from its asymptotic 



value S(t — > 00), the interblip interactions are weak and 
the WCA is expected to be a good approximation in a wide 
regime of parameters. This also implies that S(t) is not 
effective in suppressing long-blip lenghts and the NIBA 
might not be justified. Indeed, for small S s and large u) c 
(see Fig. 4a), no differences among WCA and WIBA oc- 
cur. Similarly to the Ohmic case, we show in Fig. 4b 
the parameter regime where the WCA and the NIBA are 
expected to fail. With respect to Fig. 4a, we keep here 
the same ratio w c /w p h constant, being now u c ~ A, i.e. 
the bath becoming "slow" . This case is the most difficult 
one, since the bath is very coherent and memory effects 
are to be taken into account, which requires to perform a 
very good description of the full bath dynamics. One sees 
that the NIBA completely fails to reproduce the dynam- 
ics, even reaching unphysical values. The extended-NIBA 
works better, approaching closer the QUAPI predictions. 
Nevertheless, too few correlations are taken into account, 
and it oscillates still too much with respect to the numer- 
ical plot of QUAPI. The WIBA shows discrepancies from 
the QUAPI as well, being still "too" coherent, even though 
its predictions are more accurate than the extended-NIBA. 
The WCA, despite better than WIBA in this regime, also 
lies apart from the numerical prediction of QUAPI. In 
this range of parameters, the multiphonon processes be- 
come relevant and the pcrturbative weak-coupling approx- 
imation begins to fail. This agrees with ab-initio sim- 
ulations for charge qubits interacting with piezoelectric 
phonons [11]. Hence, further analysis of the complicated 
super-Ohmic case is to be done, in order to better under- 
stand the different dynamical situations which take place 
by varying the coupling strength S s , the cutoff frequency 
u> c and the phonon frequency to p h- 
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Conclusions. — We have discussed a generalized mas- 
ter equation for the population difference P{t) of a spin- 
boson system in the whole regime of temperatures and 
couplings. This equation can be solved using standard it- 
eration schemes [17,22]. For Ohmic damping the WIBA 
is able to reproduce known results in various complemen- 
tary regimes, yielding a good, though not perfect agree- 
ment, with ab-initio QUAPI calculation in the regime of 
intermediate temperatures and damping. Hence, it over- 
comes the limits of validity of the perturbative approaches 
(NIBA, WCA) which, up to date, was possible only with 
numerical ab-initio models. For super-Ohmic damping the 
WIBA works well for large cut-off frequencies and low- 
to-modcratc temperatures. However, disagreement with 
QUAPI is found in the case of a "slow" bath. 

We mention some general contexts for the need of a 
bridging approach: i) The common experimental situa- 
tion where bath temperature or TLS asymmetry are varied 
over a wide range (WCA and NIBA are unreliable at high 
temperatures and intermediate bias, respectively), ii) Sev- 
eral TLS's interacting with a common heat bath, as e.g. in 
glasses at low-temperatures [6,7]. Due to the wide distri- 
bution of tunneling parameters and asymmetries, neither 
the WCA nor the NIBA can describe the dynamics of the 
whole ensemble consistently 

We must, however, notice that in the case of "slow" 
environments with cut-off frequency of the order of the 
tunneling frequency, our model needs some improvements, 
since neither the WIBA nor other analytical approxima- 
tion schemes are able to reproduce the correct onset of 
decoherence which in fact takes place. This situation oc- 
curs e.g. in non-adiabatic electron transfer [4] or for charge 
qubits interacting with piezo-electric phonons [11]. 

* * * 
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